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Abstract 

We compute for the first time QCD corrections to the rare decays B —> Xd,sVi>, 

K —> TTvi? and Kl —> , where f = e or in the context of a supersymmetric extension 

of the Standard Model (SM) with minimal flavour violation and new operators, in addition to 
those present in the SM. Assuming that the gluino is heavy, we consider an effective theory which 
consists of charged and neutral Higgs particles, charginos and squarks. We evaluate the QCD 
corrections to box and -penguin diagrams with top-quark, charged Higgs boson, chargino and 
squark exchanges, as well as to neutral Higgs boson penguin diagrams. We provide a compendium 
of analytic formulae for the Wilson coefficients, which are valid for arbitrary values of tan /3 (the 
ratio of the vacuum expectation values of the two Higgs fields) except for the case of the neutral 
Higgs-boson contributions. These contributions have been obtained at large tan j3, which may 
compensate for the inevitable suppression by the masses of the light leptons in decays based 
on the b —> s{d)Dl~ transition. We investigate the dependence of the various branching ratios 
on the renormalization scale p, which is the main theoretical uncertainty in the short-distance 
calculation. We find that the p dependence of the branching ratios is considerably reduced once 
the QCD corrections are taken into account. The contributions of new operators are found to 
be dominant at large tan/3 in B^^s while they are subleading in .B ^ and 

completely negligible in kaon decays. 
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1 Introduction 


The rare deeays B Bd,s —^ and Kl where I denotes a 

lepton, are very promising probes of flavour physios within the Standard Model (SM) and possible 
extensions, since they are governed essentially by short-distance interactions. (For recent reviews 
of rare K and B decays, see Refs. [M .) These decays are sensitive to the quantum structure of 
flavour dynamics and can at the same time be computed to an exceptionally high degree of precision. 
This is due to the fact that (a) the required low energy hadronic matrix elements are just the matrix 
elements of quark currents between hadron states, which can be extracted from semileptonic decays, 
including isospin-breaking effects [@]; (b) other long-distance contributions have been found to be 
negligible [|^]; and (c) the contributions of higher dimensional operators turn out to be tiny in the case 
of Kl —>• and are below 5% of the charm contribution in the case of [j|]. 

As a consequence, the scale ambiguities inherent in perturbative QCD constitute the dominant 
theoretical uncertainties present in the analysis of the above decay modes. Within the SM, these 
theoretical uncertainties have been considerably reduced through the inclusion of the next-to-leading 
order QCD corrections 

In the decay Kl on the other hand, there are sizable long-distance contributions as¬ 

sociated with a two-photon intermediate state, in addition to the short-distance interaction ds —> 

. While the absorptive contribution with real photons can be determined by means of the mea¬ 
sured Kl 77 branching ratio [p^, the numerical predictions for the dispersive part of the ampli¬ 
tude due to off-shell photons suffer from theoretical uncertainties inherent in the models for the form 
factors [|T3|]. We therefore content ourselves with the short-distance part of Kl —> . 

In the context of the minimal supersymmetric standard model (MSSM) []T^,p^], the branching 
ratios for rare 5-meson and kaon decays such as B ^ Xguu, Bg Kl K^ 

and Kl K' have been considered by many authors (see, e.g.. Refs. [|TE|-^); however, 
these decay modes have hitherto been studied without the inclusion of QCD corrections. Working 
in the framework of the SM operator basis, a compendium of the relevant formulae in the limit of 
degenerate squarks of the first two generations, and valid only for small tan (3 (the ratio of the vacuum 
expectation values of the two Higgs fields), has been presented in Ref. [0. Taking into account 
all presently available constraints on the supersymmetric parameters, the phenomenological analysis 
of PT| ] shows that the supersymmetric contributions to the above-mentioned processes can be quite 
substantial. Indeed, denoting the MSSM prediction for a given decay normalized to the SM result 
by R, and setting all the SM parameters that are unaffected by supersymmetric contributions at their 
central values, one finds the following rangesQ 


0.73 ^ R{B Xgiyp) ^ 1 . 34 , 0.68 ^ R{Bg ^ 1 - 53 , 


( 1 . 1 ) 


'Note that these ranges become larger once the SM parameters are varied. 
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0.65 ^ R{K 


TT^Uu) ^ 1.03. 


( 1 . 2 ) 


TT+i/i/) < 1.02, 0.41 ^ R{Kl 


Clearly, in view of large parametric uncertainties related to the SM parameters and the masses of 
the supersymmetric particles, precise predictions for the above branching ratios are not possible at 
present. Nevertheless, the situation may change considerably in this decade due to the improved 
determination of the SM parameters in forthcoming B- and it"-physics dedicated experiments, as 
well as through an anticipated discovery of supersymmetric particles. Besides, all the branching ratios 
considered in this paper could conceivably be measured in this decade. We think that these prospects, 
together with the theoretical cleanliness of these decays, justify a more accurate computation of the 
relevant branching ratios in the MSSM. 

In this paper, we extend existing calculations in two ways: 

(i) We compute for the first time QCD corrections to the branching ratios of B ^ Bd,s —^ 

Ti^vP, Kl —> and within supersymmetry (SUSY). 

(ii) We also include the contributions of new operators beyond those present in the SM. 

While such an enterprise would have been a formidable task ten years ago, it is a relatively straight¬ 
forward computation by means of analytic computer programs developed during the last five years. 

It is the main objective of the present work to investigate the dependence of the various branching 
ratios on the renormalization scale /r within SUSY, taking into account QCD corrections. As we 
shall see, the inclusion of these corrections will significantly reduce the unphysical renormalization 
scale dependence of the branching ratios originating in the scale dependence of the running quark and 
squark masses. 

In the quantitative analysis, we focus on the two distinctly different regions 

2 ^ tan/) ^5, 40 ^ tan/) ^ 60, (1.3) 

which we refer to as the low and high tan f3 regime, respectively. Our results for the Wilson co¬ 
efficients that will be presented below are valid for arbitrary values of tan (5, except for the neutral 
Higgs-boson contributions, which are sizable only in the high tan /) regime for the decays under study. 

The outline of this paper is as follows. In Sec. ||, we touch on the elements of the MSSM and define 
our notation. Section ^ is a compendium of the relevant effective Hamiltonians and the corresponding 
Wilson coefficients including 0{as) corrections, in the context of the MSSM. In Sec. ^ we give 
technical details of our calculation of the QCD corrections. In Sec. |^, we present a collection of the 
branching ratios in question. Section is devoted to the numerical analysis, which is based on the 
assumption of minimal flavour violation and performed in the low as well as the high tan /) regime. 
Some comments on the SUSY QCD corrections to the b quark Yukawa coupling, relevant at large 
tan (3, will be made in Sec. ^ Finally, in Sec. |^, we summarize and conclude. A compilation of the 
various loop functions is given in the Appendices. 
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2 Couplings and mixing matrices in the MSSM 


We start by introducing the relevant mass and mixing matrices in the context of SUSY with minimal 
particle content and i?-parity conservation, which we will call the minimal supersymmetric stan¬ 
dard model (MSSM). In the numerical analysis, we confine ourselves to a version of the model with 
minimal flavour violation, i.e., we assume that flavour mixing is due exclusively to the Cabibbo-Ko- 
bayashi-Maskawa (CKM) matrix. In addition, we take the down squark mass-squared matrix to be 
flavour diagonal, so that there are no neutralino contributions. 

The MSSM is usually embedded in a grand unified theory (GUT) in order to solve the problem 
with too large contributions to flavour-changing neutral currents and to further reduce the vast number 
of unknown parameters. This leads to the minimal supergravity inspired model in which one assumes 
universality of the soft terms at the scale of gauge coupling unification, say, Mgut- Renormalization 
group effects then induce flavour off-diagonal entries, e.g., in the squark mass-squared matrices at 
the electroweak scale. In our numerical analysis, however, we do not relate the soft SUSY-breaking 
parameters to some common high scale, but rather take them at the electroweak scale while discarding 
flavour off-diagonal terms in the squark mass-squared matrix. 


2.1 Chargino mass matrix 


The chargino mass matrix is given by 

M 2 y/2M\Y sin (3 
\/2Mw cos (3 n 

where M 2 and /i are the W-’mo and Higgsino mass parameters, respectively. This matrix can be cast 
in diagonal form by means of a biunitary transformation 




U*M^V^ = diag(M^^,M 


Mv. . being the chargino masses with M?^ < M? . The analytic expressions for M 


Xl,2 


X2 


Xl,2 


( 2 . 2 ) 

and the matri¬ 


ces U, V can be found in Ref. [29]. 


2.2 Squark mass matrix 

It proves convenient to work in the super-CKM basis [j^], in which the quark mass matrices are 
diagonal, and both quarks and squarks are rotated simultaneously. Denoting up-type squarks by U, 
the 6x6 mass-squared matrix in the (Ul, Ur) basis is given by 

M^^ -f M^j + M| cos 2/9(2 — I sin^ dw)"^ Mu{A^ — n cot /911) 

[Mu{AIj — jj, cot /3]1)]^ ^Or ^ + 1^1 cos 2(3 sin^ 9w^ 

(2.3) 
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where 9w is the Weinberg angle, Mfj^ ^ are the soft SUSY breaking up-type squark mass matrices, 
Mjj = diag(m„, me, rrit), Ajj contains the trilinear soft SUSY-breaking parameters and 1 represents 
the unit matrix. The matrix M? can be diagonalized by a unitary matrix such that 

= diag(m?^,m?2,... ,m?J. (2.4) 

For later use it is convenient to define the 6x3 matrices 

= (r^)e,+3. (2.5) 


2.3 Slepton mass matrices 

Similarly to the mixing matrix for the squarks, we define the matrices F^ in the charged slepton sector 
through 

= diag(m? , m|,.... (2.6) 

with the mass-squared matrix, in the super-CKM basis, 

Mf -f M|; — 2M| cos 2/3(1 — 2 sin^ 014/)]! Me{A*i — ^ia.n j3t.) 

[Me{A1 — /U tan/51)]^ -f Me — M| cos 2/5sin^ Ow'^ 

where ^ are the soft SUSY breaking charged slepton mass matrices. Me = diag(me, m^, m^) 
and Ai contains the soft SUSY breaking trilinear couplings. As before, we introduce the 6x3 mixing 
matrices 






( 2 . 8 ) 


Finally, if we assume the neutrinos to be massless, we are left with the 3x3 mixing matrix F^ in 
the sneutrino sector, which is defined by 


= diag(m? m? m? ), 


(2.9) 


where the sneutrino mass-squared matrix reads 

M? = -f cos 2/511. 


( 2 . 10 ) 
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2.4 Interactions within the MSSM 


Recalling 

= (1 =F75)/2, g = e/smOw, g] = 

M;7 = diag(m„,mc,mt), = diag(me, m^), = diag(mrf, mb), (2.11) 

the relevant interaction vertices of the MSSM with decoupled gluinos and massless neutrinos in the 
super-CKM basis can be written asQ 

^ [cot/3(-uMb/VcKM-PL<^) + tan/3('uVcKM^D-PR<^)]-f^^ + H.C., (2.12) 

yzMw 


= E 1 + V^’-Pr)' + fi'(Af‘-FL + X'r’->Pn)d] 

i=l ^ 

+ + Af»P„)Rl| + H.C., 


PiSZ = E XTimiU’.PL + V:,V„Pn + cos29„,i„)i-Z2, 




*J=1 

ig sec 9w 


{T^^T^^^)ab-^sm^9wSab 


(«: Ub)Zl 


^4 = C^uuu = --gl{u*PuT%)\ 


where 




M, 


u 


y/2Mw sin /?_ 


VcKM, 


= gaYU,2T^^V, 


M, 


D 


CKM- 


y/2Mw cos (3 ’ 

X^^ = -gVlV^, X^^ = gU,2T 


TV XIe 


A'f ■■ = -9 


\/2Mw cos (3 ’ 
Me 


U* v^L _ m ^ 

y/2Mw cos (3 


, Xf« = 0, 


(2.13) 

(2.14) 

(2.15) 

(2.16) 

(2.17) 

(2.18) 

(2.19) 

( 2 . 20 ) 


^For simplicity of presentation, we suppress the generation indices. Furthermore, we dehne the charginos (i = 1,2) 
as particles, contrary to Refs. [001111’ with = {x^ 


~-\c 
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( 2 . 21 ) 


In these equations, V"ckm is the CKM matrix and T'^ {c = 1,..., 8) are the generators of eolour gauge 
symmetry. The mixing matriees U, V, and are defined through Eqs. ([2.2D, ( p3| ) and (JZ.SP, 
respeetively, and the 3x3 sneutrino mixing matrix is defined via Eq. (|2.9|). 

The effeet of the deeoupled gluino, with mass Mg, is eontained in the funetions Ug and ay & 
whieh are given by 


= 1 - 


(/^) 

47r 


2 + 21nf4 

3 \Ml 


Cty — 1 + 


(/^) 

4:71 


l + 21n(^ 


( 2 . 22 ) 


yU being the matehing seale whieh will be diseussed in more detail at the end of See. 
As for the eouplings of the neutral Higgs bosons, we refer to [Q. 


2.5 Tree-level relations within the MSSM 


In ealeulating the neutral Higgs-boson eontributions, we shall exploit the tree-level relations 


= (2.23) 

= \{M\o +MIt [{M%> + M|)2 - 4M|Mlo cos' 2/3]'/'}, (2.24) 


sin 2a 


- sin 213 


U^o-M^J’ 


(2.25) 


where Mh and M^o are the masses of the eharged and CP-odd Higgs boson, respeetively, and M/jOpo 
and a are the masses and mixing angle in the CP-even Higgs seetor. Thus, at the tree-level, the Higgs 
seetor of the MSSM is deseribed by merely two free parameters that we ehoose to be tan (3 and Mh- 
Sinee the results that will be presented in the next seetion have been obtained at leading order in the 
eleetroweak eouplings, the tree-level relations are quite adequate for our purposes. 


3 Effective Hamiltonians and Wilson coefficients 

3.1 Notation 

Within the eontext of the SM and the two-Higgs-doublet model (2HDM), as well as the MSSM, at 
low tan (3, only a small number of operators eontribute to the deeays under study, and the Wilson 

^Equation (|2.22D can be obtained from the result of Ref. by performing the limit Mg oo. 
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coefficients can be expressed in terms of only two functions usually denoted by X and Y [|^^-|^ 
(see Appendix ^ for the SM result). 

For large values of tan (3, on the other hand, significant contributions of new operators are con¬ 
ceivable and such a convenient notation is no longer applicable. In fact, since a given Z°-penguin 
or box diagram with a certain particle exchange (IF^, x 
several operators, it is useful to proceed as follows: 


± ZT± contribute to the coefficients of 


• All Z°-penguin contributions are included in the functions and which describe the 

processes b —> qi/v and b —> (q = d, s), respectively. Likewise, the box-diagram contri¬ 

butions are expressed in terms of the functions and B^'K 

• In the case of the b ql^l~ transition, there may be sizable contributions due to neutral Higgs 
boson penguin diagrams at large tan (3, leading to the function 

• Owing to quartic squark vertices, there are additional contributions to -penguin, box- and 
penguin diagrams with neutral Higgs bosons. Indeed, when the running mass of the squarks 
in the modified minimal subtraction (MS) scheme is used, diagrams involving quartic squark 
couplings have to be included. As we shall discuss in Sec. these contributions vanish after 
renormalization when the pole mass of the squarks is used instead. Below, these contributions 
will be labelled by the index J = 4. 

• The results obtained in the B sector can be easily altered to apply to the kaon system. 

• Including QCD corrections, and denoting the SM, charged Higgs-boson, chargino and ‘quartic’ 
contributions by J = SM, H, x, 4, the functions mentioned above have the structure [| 

[cfh = \cr]f + ^[CrfX [C% = [Cfiy* + (3,1a) 

W]j = IBfl!?’ + = [Bfl!?* + (3,1b) 

[iVflj = lAlfl™ + (3,1c) 


where the subscript / indicates the various operators, which will be discussed shortly. Note 
that there are no sizable neutral Higgs-boson contributions within the SM as they are invariably 
suppressed by the light fermion masses, and hence can be safely neglected. 

‘^Recall that in our scenario (i) the gluino decouples, and (ii) there are no neutralino contributions since we assume the 
down squark mass-squared matrix to be flavour diagonal. 
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• Finally, we define the mass ratios 


X = 


K,’ 


y = 


mt 

Ml' 


z = 


Ml 

M^w' 


(3.2a) 




M\ 


9 

m~ 

2 ’ “ j^2 ’ ^>>i — -JJT ’ 


mt- 


Xi 




m 


Vfi 




Mi: 


(3.2b) 


and introduce the abbreviations 


= sin6'H^, Li = In ( -^ 1, Lu^ = \n. 

mf 




mt- 


Kr, = 




(3.2c) 


3.2 The h s{d)uvf transition 

Let us start with the effective Hamiltonian describing the b transition, where q = d,s and 

/ = e, ii,T. Within the SM, it is given byQ 

4:Gf cy 


q-fj — 


VaV-cUm.mXLbmrPLi'f), 


(3.3) 


a/ 2 27rsin^6*vi/ 

'' l=U,C,t 

where c\{mi) represent the contributions from internal quarks i = u,c, t, and rrii are the correspond 
ing masses. Using the unitarity of the CKM matrix, the effective Hamiltonian can be written as 

AGp a 


q-ff — 


a/ 2 27rsin^6'i^ ^ 

l=C,t 

AG F cy 


^ibKq[climi) - cl{mu)]{q-ft,PLb)iiyn'"PLi^f) 


y/2 271 sin^ 9w ^ 


5^ {qii^Pib) {ufYPLi^f), 


(3.4) 


i=c,t 


with c)^(O) = 0. [Note that in the second line of Eq. ( 3.4| ) we have taken the limit mu 0.] 
Exploiting the fact that 

clirnl) 


ciim) 




VM* 

* CO* cq 




0(1) {q = d,s), 


(3.5) 


the charm-quark contribution to the quark-level process b —> qi^fUj can be safely neglected, and we 
arrive at 


nL = l/bU44(mt)(g7^PL6)(i//7'^PLi//). 


y/2 2713111^9 


(3.6) 


w 


^Throughout this paper, we suppress the neuttino flavour index / carried by the Wilson coefficients. 
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The extension of this result to a non-standard operator basis is straightforward. Henceforth, we will 
omit the superscript t of the Wilson coefficient in Eq. (|3.6|). 

Assuming that the neutrinos are essentially massless, and hence purely left-handed, the effective 
Hamiltonian in the presence of new operators has the particularly simple form 

° y] [ciC>i + CHO„], (3,7) 

V2 27r sm^ Ow \ 

/=e,Ai,r 

where 

Ol = {q'yf,PLb){9f'y^PL^f), Or = (g7^PR6)(i//7'"PLZ//), (3.8) 

and the short-distance coefficients 

E {\CT]j + \BT]j}, C„= {[Cth + lBnj}- (3,9) 

J=SM,H,xA J=SM,J7,x,4 

Note that the Wilson coefficients of tensor and scalar operators such as PL^Rh){i)a^^PL,Rv) and 
{qPL,Rb){vPL,R^) vanish for massless neutrinos in the final state. 


3.2.1 Z°-penguin contributions 


The QCD corrections to the Wilson coefficients are calculated from the Feynman diagrams in Figs. 
and while those involving the quartic squark couplings are depicted in Fig. Recalling the defini¬ 
tions in Eqs. (|3^), we find 




fi^\x) + 8x-^fi^\x)Lt 


[CWhu = 0, 




L \H 


Mjj cot^ (5 
8M^ 


vPiv) + ^ 


a.. 


mfemgtan2/3 r (0). . a. 


yf2^\y) + 8y-^[yf2°\y)]Lt 
f^^\y) +8(^1+y-^y^^\y)L, 


(3.10) 

(3.11) 

(3.12) 

(3.13) 


[crh = E 


^ i^j=l a,6=1 


X 


^i\ ff\xjhyai) + ^ 


a.. 


fi^\xji,yai) + ^ [l+y, 


- ^f4\xji,yai) + 

+ |/4 \yai, yu) + 


Otg 

47r 

Otg 

47r 


fi\xji,yai) -f 4 ( 1 -f y, 


d 


d 


dya 


dyai 

c(0) 


/s {Xji,yai)Lua 


UjlU:^5ab 


yai)Lua 


v;^v,r5ab 


d 


fi^\yai, yu) + 4 (1 -f Vai^— + yu 

oya 


'dybi 


fi iqjaii yu)L{, 


(3.14) 
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V V V V 


Figure 1: Penguin diagrams contributing to the transition b —> qui? (q = d, s) at order The vertical 
and curly lines denote bosons and gluons, respectively, while the coloured particles (i.e. quarks and 
their superpartners) are represented by double lines. The diagrams for b ql^l~ may be obtained 
by replacing u I and by taking into account neutral Higgs and would-be-Goldstone bosons, in 
addition to the boson. The corresponding symmetric diagrams are not shown here. 
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b 


Q b Q 




Figure 2: Feynman diagrams with quartie eouplings eontributing to the deeay b qvu, q = d,s, 
where the vertieal line eorresponds to the boson. In the ease of the b ql^l~ transition, one 
simply replaees i/ —> /. We refrain from showing the symmetrie vertex eorreetions here. 



b Q 



Figure 3: The eontributions to the b —> quP (q = d, s) transition due to quartie squark eouplings. The 
dashed lines denote sealar quarks while the solid lines represent eharginos, sealar leptons, and the 
boson. As for the diagrams deseribing the b ql^l~ transition, we refer to the eaptions of Figs. |T] 
and 0. 
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(3.15) 


[Cr]x = [cf]x u ^v*]V ^ u*] ^ , 


[ci% = 


as KqC 


4^3M2 S 53 (Af‘'),d(Af‘).a[(fi,),t!/«(fcr)*,(l + iiJl 

iJ =1 a,...,e,g,k=l 


X ^ jifQ Uaij VdijUjiU^^dagSgi^ is^jiiyaiiydi)^j\^il^ae^gd 


( 0 ), 


+ f5^\yai, ybi, yci)(J'^'^T^^'’)bcSij6aeSbgScd + ff\yau yd, l/di)(r^^r^^'^)fec5ii(5afe(5ce<5dg (3.16) 


c(0) 


= [ci% x^«; u ^v*-v ^ u*] ^ , 

where the index q earned by the matriees is 


Q = 


1 for b ^ d, 

2 for 6 —i> s, 


(3.17) 


(3.18) 


while Pu has already been defined in Eq. ( |2.21D . We should stress that all quark and sealar quark 
masses entering the above formulae, as well as the expressions that follow, are running masses in the 
MS seheme, i.e., nib = i^biy), and so on (see See. |336| below). The loop funetions, 

/p°\ are listed in Appendix]^. 


3.2.2 Box-diagram contributions 

The relevant Feynman diagrams, which contribute at 0{as), are depicted in Figs. ^ and ^ We find 


[SflsM 


= -/f A) - 


Att 


+ Sx^fi^\x)Lt 


[^R ]sM ~ 0) 

mn = 0 , 

[B_ 




tan^ I (0) a, 

’f2(y) + ^ 


16MJM2,, 


/“(!/)+8(l+9|^)/ffe)i. 


(3.19) 

(3.20) 

(3.21) 

(3.22) 


2 2 6 , 
KnSu/ 1 


[Bfli = E E TfYip'-MPXdxfX.X'nJ /f.- m) 


2 ^ 

2 ,j=l a, 6=1 


dvr 


f8\^ji^yai,Zbi) +4(1+?/, 


d 


dya 


fb i^ji, yd, Zbi)Lu^ 


(3.23) 
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V V V V V V 


Figure 4: Box diagrams contributing to the decay b quP, q = d,s, at 0{as). The curly and 
double lines represent gluons and quarks (or squarks), respectively. The remaining internal lines 
denote gauge bosons, leptons, charged Higgs bosons, scalar leptons and charginos. The diagrams for 
b ql^l~ can be obtained by the replacements I ^ v and I z>. The symmetric diagrams are not 
shown here explicitly. 


i,j=l a,6=1 Vi f 

I r .V. I A f -I I ^ ^ .V. 1 r _ 1 \ 


“1“ /g i^ji-tyai-i^hi) H“4( 1 + tjai Zbi)L{^ 


(3.24) 


^ ij=l a,...,d=l Xi 

X [{Pu)adyM{PuUil+L^,)]f^^^ {xji,yai,ybi,Zci), 


(3.25) 


2 2 6 


a, AKaS 


ij = l a,...,d=l Xi 

X [(yPu')adydi{Pu')dbi^ P {xji-, yaii ybii ^ 


(3.26) 


where, as before, f = e, y,,T and g = 1 (2) for the 6 —d (6 —s) transition. 


3.3 The b s{d)l'^l transition 

The part of the effective Hamiltonian describing the transition b qPl~ (q = d, s) that is relevant 
for the decay Bg Pl~ has the form 


BeS — --. 2 a ^16^4[caOa + + ^sOg + c'^Og + CpOp + dpO'p], 

V 2 27 rsm dw 


(3.27) 
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where 


Oa = (qYPLb) ((7^750, (^A = (qYPRb) , 

(3.28) 

Os = mb{qPRb){ll), Og = mg{qPLb){ll), 

(3.29) 

Op = mb{qPRb)(J-f5l), Op = mq{qPLb){l'-^^l). 

(3.30) 


Again, the effective Hamiltonian [Eq. (|3.27|)] has been obtained by exploiting the Glashow-Iliopoulos- 


Maiani (GIM) mechanism [|3^], as discussed in Sec. 

As far as additional operators like and PL,Rh){la^yPL,Rl) are concerned, 

they also receive contributions within SUSY. However, as we shall see in Sec. |5.1.2| , the hadronic 
matrix elements of these vector and tensor operators vanish. Moreover, the corresponding Wilson 
coefficients do not mix with the remaining short-distance coefficients. 

The Wilson coefficients can then be written as 


J=SM,J7,x,4 


Ca = 


J=SM,J7,x,4 


(3.31a) 


J=H,xA 


c's= E 

J=H,xA 


cp= E 

J=H,xA 


c^p= E 

J=H,xA 


(3.31b) 

(3.31c) 


There are several points to be made here, (a) The SM contribution to scalar and pseudoscalar operators 
is suppressed by mirrib^q/M^, and so can be neglected, (b) The boson as a vector particle does 
not contribute to scalar and pseudoscalar operators, (c) Because of their scalar-type couplings, the 
neutral Higgs bosons do not contribute to the Wilson coefficients ca and c)^. 


3.3.1 -penguin contributions 

Evaluation of the diagrams with exchange results in 

[cih =-[Cf]j, [C%]j = -[C^Rij (J = SM,Er,x,4), 

with the functions [Cf]j and as given in Eqs. (CT)-(CT). 


(3.32) 
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3.3.2 Box-diagram contributions 


KIsm = j|/fW + S 


[^5,p]sm = [^5',P']sM = 


fil\x) + 8x-^f^^\x)Lt 


[B1]h = 0, 


[B%]h = 


mbrnqm^ tan^ P ( „/Q) , , , a* 

^n’{y) + ^ 








(3.33) 

(3.34) 

(3.35) 

(3.36) 

(3.37) 


_ m,tan^/3 / (0) a. 


a _ mztan^/3 / (0) a. 


fn\x, z) + 8x-^fi^\x, z)Lt 


fH\x, z) + 8x-^jf\x, z)Lt 


(3.38) 

(3.39) 


[Ba]x - ^ 

*j=i /=i “=i 
O^s 

Att 


" " ' (At‘'),.(.Y''‘),3 r 1 




fn’\xji,ya,Vfi) 




d 


fb {^jijyaijVfi)Li 




^ji'i /e i^ji^yai-iVfi) + 


CXg 

47r 


/9^^(a^ii,l/ai,'y/i) + 4 ( l + y, 


d 


dya 


fb yaii'^fi)Li 




Ne 


(3.40) 


[B%U = -[B%{X^^ ^ ^ X^^), 


(3.41) 


[^H = ±^EEE 


2 3 6 




mb 


i,j=l f=l a=l 


Olg 

47r 


fl2 {Xjh yah Vfi) + Ayai-^jf’ {Xji, yah Vfi)Li 

^Vai 


H 

d 


/s iXji,yahVfi) 


c(0). 




i y/Xji /g iXji,yai,Vfi) + 


tto 


fS{xji, yah Vfi) + Ayai-^f^'’’{xji, yah Vfi)Li 

^Vai 


d 


p(0) I 


(3.42) 
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(3.43) 


= ±|B!ip]j(-’f« A'‘'»i A'^*- « A"«; mt ^ m,). 

Here I denotes the external leptons (not summed over) and / labels the flavour of the intermediate 
neutrino and its sealar partner. The loop funetions fp^\ are given in Appendix p. 

Turning to the eontributions from the quartie eouplings, we obtain 

/I 2 2 3 6 ('’Writ's ( 

= E ’ .i ^" l(Pu)^ydPuUl + L..)] 


i,j=l /=1 a,h,c=l 






( 0 ), 


[B%1 = ^ ^ X^^), 


(3.44) 

(3.45) 


2 3 


[B|pl4 = T 




CK.s 4 k„Sw ^ (TD \ (-[ \ T M 

2^2^ 2^ ^- TT 2 - [{Pu)acyci{Pu)cb{l + U,)] 


All 3mb ^ ^ 

*J=1 /=1 a,b,c=l 


H 


X i f9'’i^jhyai,ybi,Vfi){X-"^^)if{X'^^)fl ± ^if[^{Xji,yai,ybi,Vfi){X-"^'^)if{X'^^)fi }>, 




[sf, p,]^ = ±[bIp]^{X^^ ^ ^ ^ m,). 


(3.46) 

(3.47) 


3.3.3 Neutral Higgs boson penguin diagrams 

As we have already mentioned, the eontributions of neutral Higgs bosons are relevant only for large 
values of tan [3. In this ease, we obtain 


tU 1 _ j ^^(0). , a.. 


mp]H = p 


4M^ 


xfs + 


dvr 


z) + 8 ^^ [ 3 ^/ 3 °^ ^)]Lt 


(3.48) 


[P^s',p']h ~ P\.^s,p]h^ 


(3.49) 




[^s,p]x ~ 3= 


Mw{Mjj - Ml 


2 'j _ _ 

i,j=l a,b=l m,n=l 


j20'Y 


X 


+ ai tan (3 + —^a 2 m^ tan^ P 
Att ^ 


Q? c 


(3.50) 


~ 3=[A^5^p]^(mq —>• rUb, X^nn Xnm)i 


(3.51) 
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where 


■pa ^ _ 

imn = ^ 

^tb^tq 


[■\/2MwVii{T^^^)naO-g “ (^!7)nn^2(r^^^)naaY] 


(3.52) 

(3.53) 


with ttg, ay defined in Eq. ( |2.22|) . The eoefficients a^’^, ai, 02 in Eq. (3.501 ) are given by 


s,P _ 

Oq — =F 


± <! f4\xij,yaj) + p 

±\ 




d 

fil\xij,yaj) + 4yai-^ f 3 °\xij, yaj)U 


U^2Vjl 


(y.c 


fSiXij, yaj) + Ayai-^fl^’^Xij, yaj)Lu^ 

^Vai 


d 


mt 




M, 


w 


ab j p(0) 


Q^o 


/s {yai,ybi) + p 


fS’hjai, 2/bi) + 4 (1 + yai^— + yu^ ) f^iyai, ybi)L. 


d 


U* V* 

j2 


d 




dybi 


Sa 

m, 


+ 


Of c 


4pa* 

imn 




fl5 {yai)^ij^ab^mm 




ai = 


y2M, 


w 


fs \yai) + 


Cts 

dvr 


fwiyai) + 4:yai^f8^\yai)Lu^ 

^Vai 


^ij^abi 


_ j2 ^(1)/ 'ire c 

^2 c\ T\ /r /l5 Kyaij^ij^ab^rrirf) 


2Mw^irr 


with 




(3.54) 

(3.55) 

(3.56) 

(3.57) 


/=i V^^Xi 

It is interesting to note that the eontribution in Eqs. ( |3.50D and ( |3.51| ) that is proportional to ai tan^ (5 
eomes from the eounterterm of the eleetroweak wave funetion renormalization. Since it is not sup¬ 
pressed by the mass of the light quarks, it gives by far the dominant contribution in the high tan (3 
regime. (Eurther details may be found in Ref. P6|].) 

Einally, the calculation of the diagrams involving quartic squark couplings yields 

2/0 2 3 6 

" 3 sr^ sr^ , tt , 

^ X ^ X ^ X ^ T^n /t^iIt \ T ■ 

Jj2aY 


47r . a,...,e,g,k=l 




^ i^PjJ^ekykji^PjJ^kgiX P tan/3 iyaii ydi)3ij5ae^gd 


+ {^f)bc 

T Mw 


^ae^gb^cdfd iyaii ybii yci) 3 “ ^ab^ce^gdfa ydi ydi 

yaj, ydj)Ui 2 Vji ± fPiXij, yaj, ydjW-^v*^ 


5ij 


^ae^dg 


{Pu)ae[l + - ff^{yej,ygj)]{Pu)gd{^t)eg5ijfr{yai,yd^ [>, 


( 0 ), 


(3.58) 
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±|]v; 


ll 1 * / \ ^ 
5,Pj4l^r 


^nm)' 


(3.59) 


Notice that there is again a tan^ (3 enhaneement in Eqs. ( |3.58| ) and ( [3.59D , due to the eounterterm 
contributions. 


3.4 The s dw transition 

The results for the s —> duV transition related to the eoupling can easily be obtained from 

the results of See. through the appropriate replaeements of flavours. However, unlike the b —> 
s{d)uu transition, the internal eharm-quark contribution to s ^ duu cannot be neglected, sinee 
\v;,vji\vivj 70, whieh partially eompensates for the suppression of the eharm-quark relative 
to the top-quark eontributions due to [of. Eq. (P3[)]. Yet, it turns out that the new-physios 

oontributions proportional to V*gV^^ are small. Aooordingly, the oharm oontribution is completely 
described by the SM for whieh next-to-leading-order oorreotions are known from [ P^[r3| ]. 

The effeotive Hamiltonian may then be written as 

W.« = “2 a Yy:.VcAO:, + V;:V„(CiOi + ChOr)], 

^J2 271 sm Ow ^ 

(3.60) 


where f = e,fi,T and 


Ol = {s'yf,PLd){uf'yf^PLUf), Or = {s'-ff,PRd){vf'y^PLVf). 

The coefficients cl and cr are given by Eqs. (|T!9D-( |3T^ if we make the replaeements 

1 


^ rudnis, Kq 


K, = 


Ss/2GFe%,Vi' 

^Ul,r 


(3.61) 


(3.62) 


Additionally, the indices of the matrices X- in Eqs. ( |3.14| )-( |3T7[ ) and ( |3.23[ )-( |3T^ should be 


changed as follows: g —> 2, 3 —^ 1. Note that the results of the previous subseotions suggest that 
the eontributions of the Wilson eoeffieient cr oan be neglected in the ease of the s ^ d transition. 
Indeed, since cr involves always the faetors msirid/M^, it is far too small to give an appreciable 
eontribution. 

As for the function in Eq. (p.60| ), it results from the next-to-leading order calculation [[TSll , 
and is given explicitly in Ref. []T2|], where numerieal values of = A and different ehoiees 


of A^ and rric may be found. Note that there is typically a 30% suppression of the eharm eontribution 
due to the QCD eorreetions. 
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3.5 The s 


dZ+Z transition 


The results for this transition can be readily obtained from those given in Sec. |^. As in the s duP 
transition, the charm contribution proportional to cannot be neglected. However, it turns out 

that the new-physics effects in the charm sector are numerically small, and thus the charm contribution 
is dominated by that of the SM. 

The effective Hamiltonian is then given by 


TdeS — 


2,G F ct 

a/ 2 27rsin^6'M/ 




A 


+ A + Ca^'a + C 5 CI 5 + Cg^'s + CpOp + c'pO'p) , 


(3.63) 


where the function Inl (the analogue of has been calculated in Refs. iBiia- The remaining 
short-distance coefficients are given in Eqs. ( [3.31D with the replacements ( |3.62| ), together with 


A. 


A _ KndK:. 

/'mn — -r r -r ! 

^td ^ts 


(3.64) 


while the operators can be obtained from Eqs. ( |3.28| )-( |330| ) through the appropriate changes of 
flavours. 


3.6 Renormalization group evolution and scale dependence 

The renormalization scale dependence of the Wilson coefficients C{fi) is governed by 


d 


dln/j, 


(7(/x) = 7^(7 (/x), 


(3.65) 


where 7 is the anomalous dimension matrix]^ Since the anomalous dimensions for all operators in 
Eqs. (pl 8 |), ( |3.28D -( ^30| ) and ( |3.61D vanish, the Wilson coefficients must be independent of /i. Other¬ 
wise the physical amplitudes would depend on the renormalization scale. 

We now demonstrate that all contributions from the -penguin diagrams, as well as from box and 
penguin diagrams with neutral Higgs bosons, are separately independent of /r when 0{as) corrections 
are taken into account. 

In the leading-order expressions of the Wilson coefficients, there are two sources of scale depen¬ 
dence. Eirstly, the running quark mass 

7/°V(2/3o) 


mg(/i) = mg(/io) 


(/^) 

cis(/^o) 


T 


w = 8, Po = n 




(3.66) 


°For a detailed discussion see, e.g., Refs. 
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where mg(/xo) is the value of the quark mass at the seale /iq and n/ is the number of active flavours. 
To first order in the above expression can be written as 


= mg(/ro) 


1 - 


Jm 


(/^O 


ln(^ 


7i°V(2/3o) 


2 dvr 

Secondly, the running squark mass 

mg(/i) = mq(/io) 

which to first order in as can be written as 

( 0 ) / 


(3.67) 


7^°) = 4 

/m 


mg(/i) = mq(/io) 


1 - 


7^ 


dvr 


ln( ^ 
/^o 


(3.68) 


(3.69) 


The /i dependence of the top-quark mass in the leading-order terms in both the SM and the 2HDM 
coefficients is cancelled by the as corrections of the Wilson coefficients that involve the derivatives 
of the leading functions [see, e.g., Eqs. ( p.lUD and ( p.l2| )]. The same applies to the /i dependence of 
the squark masses in the supersymmetric contributions [see, e.g., Eq. ( |3.14[ )]. 

As can be seen, for example, from Eq. ( |3.13|) , there is an additional /r dependence in the leading 
contributions of the new operators generated by charged Higgs-boson exchanges that is related to the 
light quark masses present in the Yukawa couplings. This dependence on /r is cancelled by the 0{as) 
corrections proportional to the leading functions. 

Einally, we observe that in the chargino sector there are additional /i dependencies in the chargino- 
squark-quark vertices, which are given in Eqs. ( |2.17[ ) and ( [2.18D . These /i dependencies arise not only 
from the /i dependence of Mu and Mu but also from the coefficients ag and ay given in Eq. ( |2.22D . 
Taking into account both /i dependencies, the effective /i dependence of the vertex is given by 


= .Y7'“(Af5) 


2^ In 

dvr 




(3.70) 


This /i dependence is cancelled by the corresponding dependence in the 0{as) corrections propor¬ 
tional to the leading functions [see, e.g., Eq. ( P.ldQ ]. We observe that the /i dependence of ag and ay 
is essential to obtain /i-independent chargino contributions. The /r dependence in Ug and ay is related 
to the fact that the effective theory does not contain gluinos, which have been integrated out at the 
scale 0{Mg). Thus, in the effective theory the gaugino and Higgsino couplings renormalize differ¬ 
ently from the ordinary gauge and Yukawa couplings. If one expresses the supersymmetric couplings 
at /i -C Mg in terms of the SM couplings such as g and the top Yukawa coupling, this difference has 
to be taken into account in order to obtain correct results. 

An additional /r dependence arises in the MSSM from the 0{as) contributions due to the quartic 
squark vertex. Note that only the /r dependence from gluonic corrections has been taken into account 
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in Eqs. ( p. 68 p and ( p.69p . The inclusion of the contribution from the quartic squark coupling results 
in a modification of the anomalous dimension 7 ^^. This additional contribution would then cancel 
the /i dependence present in the two-loop diagrams with quartic squark couplings denoted by the 
subscript J = 4 in Eqs. (3.91 ) and (3.31D . In the next section, we will give a recipe how to avoid the 
appearance of the quartic squark coupling using an on-shell prescription for the squark mass. In this 
sense Eqs. (3.68D and (3.69D can be considered as complete. 


4 Details of the calculation 


In this section we present the details of the calculation. Readers who are not interested in these tech¬ 
nical issues should proceed to Sec. ^ 

Eor processes that take place at scales much lower than Mw, models such as the SM, 2HDM and 
MSSM can be replaced by an effective theory by means of integrating out all particles heavier than 
0{Mw)- Our aim is to find the QCD corrections to the Wilson coefficients of the relevant operators 
in the effective theory. The operators of interest are given in Eqs. ( |3^ and ( 3.28D -( 330| ) for the 
quark-level transition b s{d)uu and b respectively. (The presentation given below 

applies tob ^ transitions and with obvious changes to 6 —> suV.) 

The simplest way of finding the Wilson coefficients is to require equality of one-particle-irre- 
ducible (IPI) amputated Green’s functions calculated in the full and effective theory. The former 
requires the calculation of box- and penguin-type diagrams. The different possible topologies can be 
found in Eigs. M Note that photonic penguin diagrams do not contribute to the processes under 
consideration. Within the SM, the neutral Higgs boson penguin contributions are negligibly small, 
since they are suppressed by the light lepton and quark masses and not enhanced by tan /3. All together 
we have 16 box diagrams in the SM, additional 20 box diagrams in the 2HDM and 4 box diagrams in 
the chargino sector of the MSSM. Eurther, the number of Z°-penguin diagrams amounts to 40 in the 
SM, additional 16 in the 2HDM and 20 in the chargino sector. Note that the topology shown in Eig. ^ 
is present only in the chargino sector. We have 240 diagrams containing neutral Higgs bosons in the 
2HDM and 64 neutral Higgs penguin diagrams containing charginos and squarks. While counting 
the diagrams, different quark and squark flavours, as well as different chargino generations, are not 
taken into account. Note that the up- and charm-quark contributions to the SM and 2HDM Wilson 
coefficients are included via the GIM mechanism, as explained in Sec. The sum over all squarks 
and charginos is shown explicitly in our results. 

The calculation has been performed within a covariant and gauge for the gluon and the 
W boson, respectively, which provides a useful check of our computation. The contributions of the 
penguin and box diagrams to the Wilson coefficients are separately gauge independent with respect 
to the gluon gauge, but in general gauge dependent with respect to the IE-boson gauge. Thus, the 
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results presented in the previous seetions are given in the’t Hooft-Feynman gauge for the IV boson. 
We have eheeked that the sum of penguin and box diagrams eontributing to the Wilson eoeffieients is 
independent of ^w 

We are interested in dimension six four-fermion operators, so the external momenta ean be set 
to zero from the very beginning. The masses of the light quarks ean be safely negleeted as long as 
they are not enhaneed by tan (3. That is, we negleet rritj, rrid^s and mi in the propagator, but keep them 
in the Yukawa eouplings when multiplied by a faetor of tan (3. Sinee we are interested in on-shell 
results, it is pertinent to ask whether we are allowed to keep light masses in the Yukawa eoupling 
while negleeting external momenta. In faet, higher-order terms in external momenta would enlarge 
our operator basis, but also give eontributions to Wilson eoeffieients of dimension-six operators after 
applying the equation of motion. However, these terms are not enhaneed by tan f3 as they eome solely 
from an expansion of the propagators, and therefore are negligible. 

Setting the light quark masses to zero in the propagator produees infrared (IR) divergeneies. We 
regularize the IR and ultraviolet (UV) divergeneies simultaneously m D = A — 2e dimensions. At 
one-loop level, the UV divergeneies in the penguin diagrams ean be removed by the eleetroweak 
renormalization of the wave funetion. We have ehosen an on-shell preseription, following the ap- 


proaeh of Ref. [pq]. The QCD renormalization, whieh beeomes neeessary at the two-loop level, is 
performed in the MS seheme. 

As mentioned earlier, the results for the Wilson eoeffieients presented in See. ^ are given in terms 
of the running MS quark and squark masses = mg(n) and = mg(/i), respeetively. Alterna¬ 
tively, one ean work with the pole masses, in whieh ease the following steps should be performed: 

Step 1 Remove the eontributions due to quartie squark eouplings (i.e. the eontributions with the 
index J = 4). 

Step 2 Make a shift from the MS seheme to the eorresponding pole masses, namely:|] 




g.(™r‘‘) 

Att 


-4 In 

3 


pole 


(4.1) 


mg(/i) = m 


<? 


X 


1-f 


IQ' 

(/^) 

Att 

Jm jm 

3 

(/i) - as{ml 

m 2/)2 j 

47r 


7^V(2/3o) 


(g = d,s,b), 


(4.2) 


mq{^) = m 




47r 


3 \ fx 


(4.3) 


^Note that only in the case of light quarks (d, s,b) it is necessary to resum the large logarithms. 
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where 


^( 1 ) = 
/ m 


404 40 .mo ^8 

-Uf, jJi = 102 - —Hf, 


(4.4) 


3 9 

and 7 ^^ , j3o are given in Eq. ( |3.66| ). Observe that the shift in Eq. ( p3| ) involves only the gluonie 
eorrections, sinee the eontributions due to quartic squark couplings have already been considered in 
step 1. In this context, we would like to remark that the absence of the ‘quartic’ contributions in an 
on-shell scheme is related to the renormalization of the squark mass and mixing angle. (Eor details, 
we refer the reader to Ref. P7|].) 

After the proper renormalization, the left-over divergencies in the full and effective theory must be 
of infrared character, and they must cancel each other in the matching procedure. Note that in order 
to obtain this cancellation, the treatment of light masses must be identical in the full and the effective 
theory. As the light masses are zero, all unrenormalized loop diagrams in the effective theory vanish 
because of a cancellation of IR and UV divergencies, which considerably simplifies our calculation. 
The UV counterterms in the effective theory must therefore exactly reproduce the IR divergencies of 
the full theory, which provides a further check of the calculation. We wish to emphasize that owing 
to the dimensional regularization of the IR divergencies, the matching procedure must be performed 
in D dimensions [PS|]. 

In intermediate steps of the calculation structures like 




(4.5) 


occur with Pa,b being either Pl or Pr. They cannot be reduced using D dimensional Dirac algebra, 
due to the appearance of the matrix 75 . Only after the matching all divergencies cancel and the limit 
D —> 4 can be taken. Consequently, evanescent operators must be introduced in the effective theory. 
These operators appear in box diagrams in the SM and in the MSSM. Eor example, in the SM we 


define the evanescent operator for b sPl as follows : 

Of = (s7„,7„,7,3PL&)(r7“^7“^7“^^L/) - A{s'y^PLb){h^PLl). (4.6) 

(Eurther details can be found in [jlT]].) In the chargino sector of the MSSM, operators with a different 
spinor ordering show up in box diagrams: 

OsAB = {sPAl)(iPBb), (4.7a) 

OvAB = {s-fM){h"PBh), (4.7b) 

Otaa = (4.7c) 


We are not allowed to project these operators onto the operators given in Eqs. ( |3.28[ )-( |J3(]t ). Eor such 
a projection we would have to apply Eierz identities which cannot be continued to D dimensions. Eor 
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this reason, we have to define the following so-ealled ‘Fierz-vanishing’ evaneseent operators 


'^SLL 

= OsLL 

+ -jOsLL — 

-^Otll, 

(4.8a) 

/AE 

'^SLR 

= OsLR 

+ -^OvRL, 


(4.8b) 

'^VLL 

= OvLL 

— OvLL, 


(4.8e) 

qE 

^VLR 

= OvLR 

+ ‘^OsRL, 


(4.8d) 

/AE 

'^TLL 

= Otll 

- QOsll - 

^Otll, 

(4.8e) 


as well as operators whieh ean be obtained by an interehange of Pl ^ Pr- The operators without 
tilde are identieal to the ones given in Eqs. (Q), but with exehanged I and b spinors. Moreover 
o'fj.u = Due to a finite mixing into the physieal operators the ‘Fierz-vanishing’ evanes¬ 

eent operators eontribute at next-to-leading order. A similar situation was deseribed in [ [TTI ] . From 
Fq. ([4.8eD it is evident that we eannot negleet the eontributions to the tensor operator Oraa as it af- 
feets the Wilson eoeffieients of the sealar operators. Furthermore the following evaneseent operators 
are neeessary at intermediate steps: 


Of = (s 7 „i 7„2 7„3 Pl/)(/ 7“17“27“3P^6) - IQOvll, 
Of = {s^a,'ya,la,PLl){h"^l^^'y"^PRb) - AOvLR. 


(4.9) 


5 Branching ratios 

5.1 Rare B decays 

5 . 1.1 B Xd^si^v 

The deeays B —> XqUV (q = d, s) are the eleanest theoretieally in the field of rare B deeays. Sinee 
the neutrinos eseape deteetion, these deeays are probed by requiring very large missing energy, Eg = 
{Er — Exfj/rrib (for a full diseussion see, e.g.. Ref. [jT^]). 

Using the effeetive Hamiltonian given in Fq. (p.7p, we obtain the differential deeay rate 

dB {B ^ Xquu) _ | P (P —Xcez/e)K(0 )-,3 - 211 / 2 '^ 

sin^ Ow Y 

X |(|clP + |crjn[(l - Eq){4:Eq - 1) +mj(l - 3Eq)] + 6mq{l - 2Eq - mJ)Re (c^c^) |, 

(5.1) 
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where the sum runs over the flavour of the three neutrinos, rhi = rrii/mb, and the various input 
parameters are listed in Table The expression for the missing energy speetrum is equivalent to the 
result presented in Ref. [p^. In writing Eq. (|5.1D, we have negleeted non-perturbative eorreetions of 
0{l/ml) rP3|| and 0(l/m^) [p4],p31] whieh have been found to be small over most of the Dalitz plotj^ 
Furthermore, negleeting non-perturbative eorreetions, we have used the inelusive semileptonie deeay 
rate 

r(S ^ X,ePe) = (5.2) 

in order to remove the uneertainties due to an overall factor of mf. The functions /(rhc) and K{rhc) 
represent the phase-space and the one-loop QCD corrections, respectively 

/(rhc) = 1 — 8m^ -f- 8m® — m® — 24m^ Inrhc, (5.3) 


nirhc) = 1 -f 


as{mb) Ao{mc) 


TT fifric) ’ 

where Aq can be found in Ref. pT}]. Expanding K{mc) in Eq. ( ^T^ around m-c = 0.3 results in 

as{mb) 


(5.4) 


K{7hc) — 1 


71 


1.670 + 2.027(0.3 - me) + 2.152(0.3 - me)^ 


(5.5) 


which is accurate to better than 1%, and hence perfectly adequate for our purposes. Further, 


k(0) = 1 + 


as{mb) 


TT 


25 2 , 

- 71 

6 3 


= 0.830 


(5.6) 


represents the QCD correction to the matrix element of 6 —> quu due to virtual and bremsstrahlung 
contributions [Q. 

Integration of Eq. ( |5.I| ) over (1 — m^)/2 ^ Eg ^ (1 — mg) then yields the branching ratio 

\VbbV*\‘^B{B^X,eu,)K{0) 


XgUP) = 


ATT^sin^ 9w \Vcb\‘^ /(mc)K(mc) 


X 




CLC*R)mgf{mg) 


(5.7) 


where 


/(mg) = 1-1- Om-g — Om^ — m® -f 12mg(l -f m^) Inm^, (5.8) 

and /(mg) is given in Eq. (^3l). 

^Unlike the b —> qAl~ transition, there is no virtual photon contribution in 6 ^ qvv, so that non-perturbative 
corrections due to charm quarks are further suppressed [Q]. 
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Table 1: Input parameters used in our numerieal analysis. The isospin-breaking eorreetions tk and 
the eharm eontributions Xnl = (2-^nl + ^nl have been taken from Refs. [0 and [Q, 

respeetively. Note that rric is the running charm quark mass in the MS scheme normalized at rric. 


For the CKM matrix, we use the standard parametrization [pd|], with four independent parameters 
-512, Sia, S 23 , S. As for the remaining parameters, we utilize the values compiled by the Particle Data 
Group [P^]. 


Quantity 

Value 

sin^ 6w 

0.23 

ots{Mz) 

0.118 

Sl2 

0.222 

Sl3 

3.49 X 10-3 

S23 

0.041 

5 

57° 

\Vtd\ 

7.81 X 10-3 

IBs 

230 MeV 

pole 

ml 

4.8 GeV 

pole 

ml 

174.3 GeV 

rric 

1.3 GeV 

a 

1/129 

ri 

1.17 X 10-^ 

r2 

0.24 

rs 

13.17 

rK+ 

0.901 


0.944 

Anl 

9.78 X 10-^ 

^NL 

3.03 X 10-^ 

BiB^XceUc) 

10.58% 


Table 2: Numerical values for the running quark masses rrii = employed in our analysis. 


Scale 

mt [GeV] 

ms [MeV] 

mb [GeV] 

pole 

/i = mf 

166 

61 

2.9 

pole 

ir = ml 

— 

90 

4.4 

/i = 2 GeV 

— 

110 

— 
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Thus far, no attempt has been made to seareh for the inelusive b duu deeay, and so we eoneen- 
trate on the b ^ s transition. The best upper limit has been set by the ALEPH Collaboration pS|]f| 


(5 ^ < 6.4 X 10“^ (90%C.L.). 

Using the numerieal values listed in Tables [I| and ||, together with Eq. (|5^ , we obtain 




XsVV) = 5.39 X 10" 


X 


6 

■ 0.53 ■ 

■ 0.88 ■ 

\B{B^X,ePe)] 




_K{rhc). 

10.58% 

iKfel 


r+icRp) 

Re [ctcX] 

1 0.08, 


L |cLr + |cKkj 


whieh will be used in the subsequent analysis. We note in passing that, sinee 

Re (c^c|j) 


1 

- < 


IclT + \cr 


1 

2 ^ 2 ’ 


(5.9) 


(5.10) 


(5.11) 


the last term in square braekets in Eq. ( |5.10| ) deviates from unity by at most 4%. The SM result is 
obtained by summing over the neutrino flavours /, and taking the limit cr —0 while retaining only 
the SM eontribution in c^. In this ease, cl = -Asm, with Xsm given in Eq. ( |A.ip of the Appendix. 


5.1.2 Bd,s 


The deeays Bq are after B —> Xquu the theoretieally eleanest deeays in the field of rare B 

deeays. In faet, like in the deeay B —> X^^gUp, the eharm eontributions are eompletely negligible. 
These proeesses, whieh are dominated by Z°-penguin and box diagrams, have been studied by a 
number of authors [|24|-|2^,|5T|,^ in extensions of the SM, but without QCD eorreetions. 


Eet us start by eonsidering the matrix element for the deeay Bq l^l in the presenee of the 
operators defined in Eqs. ( |3.28D -( ^3D| ), whieh has the general form[^ 


M = -tfs, ^ V,,V:q[Fsll + Fph,l + Fyp^h.l + FApni.l^ll (5-12) 

V 27r sm 6w 

where the Fj’s are Eorentz-invariant form faetors, is the four-momentum of the initial B meson, 
and fsq is the eorresponding deeay eonstant defined via the axial veetor eurrent matrix element 


{^\(lliilbb\Bq{p)) =ip^,fB,, (5.13) 

similar upper limit has been derived by the CLEO Collaboration for the branching fraction B {B^ K^vV) < 
2.4 X 10“"* at 90% C.L. We do not address the issue of exclusive decays here (see, e.g., Refs. @11). 

'^Observe that there is no tensor-type interaction as {Q\qaqvb\Bq(p)) = 0. In fact, it is not possible to construct a 
combination made up of that is antisymmetric with respect to the index interchange p ^ v. 
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while the matrix element of the veetor eurrent vanishes. It may also be noted that the form factor 
does not contribute to the decay Bq - 
find for the remaining matrix element 


V 


l^l since l^l = 0. Employing the equation of motion, we 


{^\qibb\Bq{p)) = -ife^ 


Ml 


(5.14) 


rrih + mq 

Squaring the matrix element and summing over the final lepton spins, the branching ratio can be 


written in a compact form pbQ: 

B{Bq 


GWMb fl 

i+r) = 


IbTT^ sin^ 9w 




dmf r 

MIX 


1 - 


4m? 


iFsl"^ + \Fp + 2 miFAp 


(5.15) 


Here, tb is the lifetime of the Bq meson and 


Fs = 


Cs - Cgniq 
1 + m-o 


Fp = -Ml 
2 


Cp — CpTUq 

1 + m 


g J 


Fa = -(ca-c'a) 


(5.16) 


(remembering that ?hg = mq/mb). At present, the best upper limit on the above decay modes comes 
from the Collider Detector at Fermilab (CDF) and has been derived for the b ^ s transition : 

B {Bs /i+/i“) < 2.6 X 10"® (95% C.L.), (5.17) 

and so we focus on the Bg F' decay. 


Recalling the scalar and pseudoscalar Wilson coefficients in Sec. |^, it turns out that the contri¬ 
butions of the operators Os,p and O'g p are of comparable size. Thus, the Wilson coefficients c'g p in 
Eq. ( |5.16| ) can be neglected since m^ -C 1 for g = d, s. 

Introducing the dimensionless Wilson coefficients 


Cs = Mb^cs, Cp = Mp^cp, 

the branching fraction is given by 

B{Bs /iV) = 2.32 X 10 


(5.18) 


-6 


^Bs 

IBs 

2 

r iEi.i 1 

1.5 ps 

230 MeV 


_0.040_ 


X [0.998|csP + \cp + 0.039(ca - 


(5.19) 


where we have set m^ = 0 in Eq. ( |5.16| ) and used the input parameters shown in Tables |I] and |^. The 
SM result for the branching fraction may be obtained from Eq. (|5.19D by setting cg = cp = 0(4 = 0 
and ca = — Esm, with Esm defined in Eq. ( |A.2D of the Appendix. 
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5.2 Rare K decays 


In this section we adopt a somewhat different notation from that of Refs, 
high powers of \ Vus\- That is, we define the ratios 


in order to avoid 


_ a^B {K^ 7r°e+t/e) tk^ _ _ B + ^ 

^ sin^ I Ks P tk+' ^ ^ B {K+^ 

with their numerical values summarized in Table [I]. 


(5.20) 


5.2.1 K+ 




Using the effective Hamiltonian in Eq. ( p.60D , it is straightforward to find the branching ratio for this 
decay. Since the matrix element of the operator Or equals the known matrix element of Or, we 


readily obtain the branching ratio from the formula given in Refs. 0121,p^: 

B {K^ = 2rir2rK+ ^{[Im AfX(a:)]^ + [RcAcXnl + Re AiX(a;)]^}, 

/ 


(5.21) 


where the sum is over three neutrino species, rx+ represents an isospin correction that one encounters 
when relating vr+i/z/ to vr^e+i/e 0, and Aj = V*gV^^. The function 




denotes the charm contributions discussed in Sec. |3)4| , and 

X = Cl + cr 


(5.22) 


(5.23) 


replaces the SM function Xsm given in Appendix ^ The numerical values of ri 2 , rx+, and Xnl are 
given in Table [I|. 

As far as the current experimental situation is concerned, the first clean X+ —vr+i/z/ event 
was found by the E787 Collaboration []5^]. Recently, further evidence for this decay mode has been 
reported in Ref. [ |55| ] . The updated branching ratio 

B {K+ 71+uu) = {l.57+lil) X 10“^°, (5.24) 

being roughly by a factor of two higher than the SM expectation, provides already a non-trivial lower 
bound on Vtd when interpreted within the SM framework [ p5| , |5^ . 
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5.2.2 Ki 




The decay iTi —is a short-distance dominated process that is largely governed by CP-violating 
contributions and thus is sensitive only to the imaginary parts of the CKM couplings.As 

the coupling V*gV^^ is real to an excellent approximation, the internal charm contributions can be 
completely neglected. 

The branching ratio may be obtained from the usual SM expression []T^[r3|] by 

replacing Asm through X as defined in Eq. ( p.23| ). Summing over neutrino flavours, we have 


B {Kl TT^z/i/) = 2rirK^ } jImAtA(a;)]^, (5.25) 

/ 

where the isospin correction factor is given in Table [^. The best current upper limit has been set 


by the KTeV Collaboration [ pSQ : 

B {Kl ttVp) < 5.9 x 10"^ (90% C.L.). 


(5.26) 


5.2.3 Kl —^ 

As mentioned at the outset of the paper, the decay Kl iKjji~ suffers from theoretical uncer¬ 
tainties due to the long-distance dispersive contribution 0], and we therefore concentrate on the 
short-distance effects. In the part oc only the SM operator Oa has to be kept, as in the de¬ 

cays K^ Ti~^uP and Kl tv^uP, since the Wilson coefficients of the remaining operators are 
completely negligible. Likewise, in the part proportional to the CKM elements non-negligible 

contributions arise only from Oa- 

Recalling A* = and Eq. ( |5.20| ), the branching ratio may be written as 


B {Kl )sD = 4rir3[Re AcTnl + ReAiE(a;)]^ 


(5.27) 


where Inl represents the charm contribution obtained in Ref. [ [T^ and Y{x) = —ca- The numerical 
values of r^, as well as the remaining parameters, are listed in Tables and ^ 


6 Numerical analysis 

In the subsequent analysis we will adopt the following procedure: 

• We restrict our attention to the low and high tan f3 regime, as defined in Eq. (|1.3[). Eor the 
decays with a dilepton in the final state, we focus on the mode. 


' * Strictly speaking, this does not pertain to models with lepton flavour violation, in which the CP-conserving amplitude 
can dominate 
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• In the case of B decays, we will study the /x dependence of the various branching fractions only 
for those decay modes that are mediated by the b ^ s transition, which essentially involves 
the CKM element Vts- It is important to emphasize that the presence of new-physics contri¬ 
butions may not only affect the decay modes under study but also B^-B^ mixing and the CP 
violation parameter ex, and consequently the extraction of the CKM elements Vtd and Vts- In 
this case, the standard analysis of the unitarity triangle may lead to false results. In fact, while 
to a good approximation \Vts\ ~ \Vcb\, independent of new-physics effects, the value of Vtd 
determined using the SM formulae might differ from that obtained in the context of SUSY (see, 
e.g.,Refs. [|2T|,|3|,|5l]). 

Since we are mainly interested in the /i dependence of the various branching fractions, rather 
than on their actual values, we fix | Vtd\ to the SM value given in Table [^. Note that this treatment 
is different from the analysis of Ref. [^Tj], where the new-physics effects on Vtd have also been 
taken into account. 


• The light quark masses, rUg and mi,, appearing in the Wilson coefficients are determined at the 
high-energy scale, but otherwise are evaluated at the low-energy scale (cf. Table |^). Since the 
contributions proportional to the down-quark mass are negligibly small, we may take = 0. 

• Our calculation is based on the assumption of minimal flavour violation, as outlined in Sec. ^ 
Furthermore, we assume that there are no new CP-violating phases in addition to the single 
phase residing in the CKM matrix. 


Since we ignore flavour-mixing effects among squarks, the matrix in Eq. ( |2.3| ) decomposes into 
three 2x2 matrices. A noticeable feature is that the LR terms are proportional to the masses 
of the up-type quarks. Hence, large mixing can occur in the scalar top quark sector, leading to 
a mass eigenstate, say, ii, possibly much lighter than the remaining squarks. We therefore keep 
LR mixing only in the stop sector, where the mass matrix is given by 


M| = 








mt{At — /X cot j3) 


+ m'f + I Af| cos 2(3 sin^ 


Uw 


( 6 . 1 ) 


where ^ are the soft SUSY breaking scalar masses and At is the trilinear coupling. In this 
framework, the mixing matrices P^^ and P^^ [Eq. (|2.5[)] take the simple form 



(l 0 

0 

0 0 

0 \ 


/o 0 

0 

1 0 

0 \ 

(P^-)T = 

0 1 

0 

0 0 

0 


0 0 

0 

0 1 

0 


\0 0 

cos 6^ 0 0 

- sin Oi / 


\0 0 

sin 6i 0 0 

cos 6i y 
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The physical mass eigenstates are then given by 


ti = cos OitL + sin OitR, t 2 = — sin Oiti + cos OiIr, 


(6.3) 


with the mixing angle {—'k/2 ^ 9i ^ 7r/2) 


sin 26 1 = 


2mt{At- cot f3) ""in' ' 6 


9 2 

m~ — m- 

C 2 


cos 26 *£ = 


(m? — m? ) +3M| cos 2/9(3 — 8 sin^ 6'vi/) 


2 2 

ml — m| 

Cl C2 


, (6.4) 


2 being the stop masses with . (The remaining up-type squark masses are taken 

to be equal.) 


• For simplicity, we assume that the scalar partners of the leptons are degenerate in mass. 


• For the results presented below we take into account the following lower bounds on the SUSY 
particle masses 


m 


il,bl 


> 


100 GeV, 


^q^h,bi ~ 


260 GeV, m,-- > 100 GeV, 


> 100 GeV. 


(6.5) 


As far as the lightest neutral Higgs boson, h^, is concerned, we must ensure that M^o ^ 
113.5 GeV [^1^], taking into account radiative corrections to the tree-level mass 

defined in Eq. ( |2.24| ). 


Further restrictions on the SUSY parameter space are imposed by electroweak precision data 
such as the p parameter []42|,|5D|,|55|]. Also, we take into account the constraint on the trilinear 
coupling. At, arising from the requirement of the absence of colour and charge breaking minima 

B- 


• We require the various SUSY contributions to be consistent with the measured inclusive h ^ s'y 
branching fraction [^. To be specific, we will allow the range of 2.0 x 10“"^ to 4.5 x 10“"^ 
for the branching ratio B {B ^ Besides, we take into account constraints on the Wilson 
coefficients arising from other rare exclusive B decays such as Bg pApL~ [Eq. (|5.17D]. 


To determine the impact of the QCD corrections on the various branching ratios, we examine their 
p, dependence for given points in the SUSY parameter space. For definiteness, we have chosen the 
following SUSY parameter sets: 


tan (3 = 3: 


my = 200 GeV, my = 800 GeV, 6i = —70°, m^-^ = 100 GeV, 
Mh = 300 GeV, p = —300 GeV, M 2 = 800 GeV, Mg = 1 TeV, 


( 6 . 6 ) 
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taii/S = 40: 


(6.7) 


mi^ = 120 GeV, mi^ = 500 GeV, 9i = —70°, = 100 GeV, 

Mh = 250 GeV, /x = -350 GeV, Ms = 800 GeV, Mg = 1 TeV, 

o I = 500 GeV, mr = 650 GeV, 6f = —44°, mj - = 100 GeV, 

tan/3 = 40: ’ (6.8) 

[Mh = 200 GeV, /x = -600 GeV, Ms = 800 GeV, Mg = 1 TeV. 

Note that in the case of 6^ = —70° the lighter scalar top quark is predominantly right-handed 
[cf. Eq. (|6.3D] while the choice = —44° corresponds to a scenario with almost maximal mixing, 
i.e. I sin(26*£)| 1. 

In Figs. we summarize our results for the /x dependence of the branching ratios of -B —X^z/P, 
Bs —> —> TT+z/P, Kl —»• 7r°z/P and Kl —> Overall we see that the dependence 

of the various branching fractions on the renormalization scale is considerably reduced once QCD 
corrections are taken into account. In fact, the /i dependence inherent in the leading-order predictions, 
typically 10-20%, is reduced to a few per cent once 0{as) corrections are taken into account. 

Figures ^ and show the /x dependence of the B —> X^z/P branching ratio for the low and high 
tan (3 regime, respectively. Notice that the SUSY contributions interfere constructively or destruc¬ 
tively with those of the SM, depending on our choice of SUSY parameters. Another noticeable 
feature is that the value of /x at which leading order (FO) and next-to-leading order (NFO) correc¬ 
tions coincide depends on whether one uses zt = 1 or zx according to Eqs. (|53|) and ( |5.6|) in the FO 
expression for the B X^z/P branching ratio. 

Turning to the Bs —^ /i’'“/x“ decay, we show in Figs. 0 and ^ the branching fraction for the low 
and high tan (3 regime, respectively. While in the former region the branching ratio in the SM is only 
mildly affected by the SUSY contributions (see Fig. in the high tan f3 regime the supersymmetric 
effects can be enormous. This can be seen from Fig. || where we have plotted the branching ratio 
for the case of a predominantly right-handed light stop quark [Fig. ||(a)] and for the scenario of 
maximal mixing in the scalar top-quark sector [Fig. ^b)]. Note that the value of /x at which the 
leading and next-to-leading order results for the branching ratio coincide depends on the choice of 
input parameters. 

We should mention that the large new-physics effects in Bg —*> /x’''/x“ in the high tan f3 regime, 
compared to those in B ^ X^z/P, are due to the fact that in the former decay mode the leading 
contribution scales roughly as ~ tan^ /3 while in the latter decay it behaves as ~ mi,ms tan^ (3. 

In Fig. 1^ we have plotted the branching ratios for the decays /x’'“/x“, X+ —vr+z/P and 

Kl —> 7r°z/P in the low tan f3 regime. For large values of tan (3 the SUSY contributions are negligibly 
small. As far as new operators are concerned, they do not give any sizable contributions, due to the 
suppression of the light quark masses. 
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Figure 5: Predictions of the /r dependence of the B branching ratio obtained (a) using n 

as given in Eqs. ( |53| ) and ( |5.6D in the leading order (LO) as well as the next-to-leading order (NLO) 
branching ratio [Eq. (|5.7[)], and (b) taking k = 1 in the LO expression of the branching ratio. The 
solid (dashed) curves represent the SUSY results with (without) QCD corrections. We have chosen 
tan (3 = 3, together with the parameter set given in Eq. (|6.6|). Eor comparison, we also show the SM 
prediction with (dash-dotted curve) and without (dotted curve) QCD corrections. 



Eigure 6: Branching ratio of the B Xsi^u decay vs the renormalization scale /r at large tan (3, as 
defined in Eq. (|0|). We have used k according to Eqs. 


and (p.6| ) in both the LO as well as the 
NLO expression for the branching ratio [Eq. (|5.7D]. The legends are the same as in Eig. 
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Figure 7: The /r dependenee of the Bg branehing ratio for (a) the SM, and (b) SUSY in the 

low tan/3 regime, as defined in Eq. (|^). The solid and dashed eurves denote the predietions with 
and without QCD corrections, respectively. 




Figure 8: The /r dependence of the Bg branching ratio within SUSY at large tan/3. The 

solid and dashed curves denote the SUSY prediction with and without QCD corrections, respectively, 
(a) For the case of a predominantly right-handed light stop quark, using the SUSY input parameters 
given in Eq. (|6.7D. (b) For the case of almost maximal mixing in the scalar top quark sector according 
to the parameter set given in Eq. (|6.8D . Note the order-of-magnitude enhancement of the branching 
ratio, compared to the SM and low tan (3 SUSY predictions in Eig. ^ 
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Figure 9: The Kl —> and Kl tv^uu branehing ratios, as a funetion of 

the renormalization seale fi. The dash-dotted (dotted) eurves eorrespond to the SM predietion with 
(without) QCD eorreetions while the solid (dashed) eurves denote the SUSY predietion with (with¬ 
out) QCD eorreetions. We have ehosen tan {3 = 3 along with the SUSY parameter set displayed in 
Eq. d^. 


36 












7 Corrections to the 6-quark mass at large tan /3 


Before summarizing, a few remarks are in order regarding eorreetions to the b quark Yukawa eoupling, 
whieh may be important in the high tan (3 regime As a matter of faet, sizable eontributions to 

the down-type quarks, and henee to the CKM elements, may oeeur. Reeently, it has been pointed out 
in Refs. [|2^,|^,|^,^] that the eorreetions to the b quark Yukawa eoupling, whieh are eomplementary 
to those presented in the preeeding seetions, ean be substantial in rare B deeays like Bg 
and b —> 57 . Sinee our main emphasis has been on the /r dependenee of the various branehing ratios, 
rather than on their preeise values, this issue will not be pursued here. 


8 Summary aud couclusious 

We have earried out, for the first time, a ealeulation of QCD eorreetions to the branehing ratios of the 
deeays B Bd,s —^ vr+i/z/, Kl and Kl within SUSY. 

Our results are applieable to a version of the MSSM in whieh (a) the gluinos deeouple; (b) flavour 
violation is governed exelusively by the CKM matrix; and (e) the neutralinos do not eontribute. 

Our main results ean be summarized as follows: 

(i) We have provided a eompendium of branehing ratios and Wilson eoeffieients in the presenee 
of supersymmetry, ineluding all relevant dimension-six operators. Our results are valid for ar¬ 
bitrary values of tan (3 exeept for the neutral Higgs-boson eontributions, whieh have been ob¬ 
tained in the high tan /? regime. 

(ii) The inelusion of QCD eorreetions within the SM and SUSY leads to a signifieant reduetion 
of the unphysieal renormalization seale dependenee, whieh is related to the running quark and 
squark masses, and whieh is unavoidably present in the existing leading-order formulae. While 
in the leading-order expression the seale uneertainty is typieally 10 - 20 % in the branehing ratio, 
it is redueed to a few per eent after 0{as) eorreetions are taken into aeeount. 

(iii) For the set of SUSY parameters eonsidered in this paper, it is possible to find a value of the 
renormalization seale /i for whieh the 0{as) QCD eorreetions to the branehing fraetions are 
below, say, 2%, so that they ean be negleeted. In this ease, one ean estimate these QCD- 
eorreeted branehing ratios merely by using the leading-order result with the proper ehoiee of /i. 
Its aetual value depends, of eourse, on the proeess eonsidered, as well as on the speeifie ehoiee 
of the SUSY parameters (see our diseussion in See. |^). 

(iv) As far as new operators are eoneerned, our analysis implies that in the low tan/3 regime 
2 ^ tan /3 ^ 5, only the SM operators are relevant, with their Wilson eoeffieients modified by 
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the presence of non-SM contributions. By contrast, in the large tan (3 region 40 ^ tan (3 ^ 60, 
the effects of new operators can lead to an order-of-magnitude enhancement in the decays 
Bd,s ^ in accordance with previous studies and also provide substantial con¬ 
tributions to AMb, As for the decays B —> the supersymmetric corrections 

due to new operators, although enhanced by powers of tan (3, are generally smaller since they 
are suppressed by the light quark masses 5 . Finally, the corresponding corrections in kaon 
decays are completely negligible. 
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A Standard model notation 


In this Appendix we give the correspondence of our results with the notation commonly used in 
the literature within the context of the SM. In this case, the decays Kl B Xguu and 

13q —> are described by the loop functions Xsm and Ism in the top-quark sector, which are 

defined as [jl], ^ 


0.0 


Xsuix) = Xo{x) + = C{x) - AB{x, 1 / 2 ), 


Of c 


>SM(a:) = Yo{x) + = C{x) - B{x, -1/2), 


(A.l) 

(A.2) 


where C and B correspond to Z°-penguin and box-type contributions, respectively. In terms of the 
loop functions given in Appendix we find 


1 


C{x) = fl°\x) + ^ 


Offi 


dvr 


fi^\x) + Sx^f^'^’ix) In ( ^ 


d 


( 0 ), 


dx 




rrit 




Os 

dvr 


fG^\x) + 8 x^/ 2 ® (x) In ^ ^ 


mi 




a.. 


fil\x) + Sx^f^'^’ix) In ( 


d 


io)i 


dx 


/i" 


mt 


(A.3) 

(A.4) 

(A.5) 


The expressions in Eqs. ( |A.3| )-( |A3[ ) agree with the results found in Refs. [ ]rT| , [T^ ]. Note that in these 
papers the explicit /i dependence is given in terms of ln(/i^/M^). 
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B Auxiliary functions 

Defining the dilogarithm Li 2 by 

Li2(^) = - 

Jo t 

the loop functions fp^\ fp}^ appearing in the formulae of Sec. ^have the following form: 




= 


a:(6 — x) x{2 + 3a;) 


+ 


2(a;-l) 2(a;-l)2 


Inx, 


X 


+ 


X 


X — 1 (a; — !)■ 


In a:, 


( 0 ). . a:lna; y\ny 

/a {x,y) = 7-7T7-r + 


{x-l){x-y) {y-l){y-x)' 


f4\x,y) = 


x‘^ In a; 


+ 


\ny 


fi^\x,y,z) = 
f6^\x,y,z) = 


{x-l){x-y) {y-l){y-x)' 

x^ Inx 


(x — l)(x — y)(x — z) 
xlnx 


(x — l)(x — y)(x — z) 


+ (x ^ y) + (x <-y z), 
+ (x y) + (x z), 


fi^\x,y) = 


X Inx 


+ 


x\ny 


{x-l){x-y) {y-l){y-x)' 


_ 3;lnx 
X — l’ 


= 


f^°\w,x,y,z) = 
fw{w,x,y,z) = 


Intu 


{w — l){w — x){w — y){w — z) 
w\nw 

{w — l)(tu — x){w — y){w — z) 


+ {w x) + {w ^ y) + {w z), 
{w ^ x) + {w y) + {w ^ z), 


fihx^y) = + / V 

(x -y) [y- x) 

.(i), , 4x(29 + 7x + 4x^) 4x(23 + 14x + 3x^) , 4x(4 + x^).^, / 1 , 
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